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. Abstract 

The dilute Heisenberg ferromagnet is studied taking into account fluctuations of magnetization caused by disor- 
der. A self-consistent system of equations for magnetization and its mean quadratic fluctuations is derived within 
the configurationally averaged two-time temperature Green's function method. This system of equations is analiscd 
at low concentration of non-magnetic impurities. Mean relative quadratic fluctuations of magnetization are revealed 
to be proportional to the square of concentration of impurities. 
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Introduction 



The method of configurationally averaged Green's function developed by Kaneyoshi |l|, |2j to study a dilute ferromagnet 
was used successfully by several authors [|[ [|, ||, |(| to describe disordered magnetic systems. One of the most frequently 
used approximations within the frames of this method is neglecting the fluctuations of magnetization caused by disorder 
[Q H . Note that at zero temperature the magnetization for Heisenberg model tends to saturation and fluctuations are 
suppressed. 

In paper [Q it was shown that taking into account the structural fluctuations of magnetization is necessary in order 
to explain the anomalous behaviour of spin-wave stiffness constant in amorphous ferromagnets. It was assumed that 
the structural fluctuations of magnetization at each site are statistically independent and are governed by a Gauss 
distribution function. In order to obtain quantitative results one needs to take into account the connection of the 
structural fluctuations of magnetization with the fluctuations of structure. In our previous work || this relation was 
obtained within the framework of mean field approximation. It allowed one to analyse quantitatively the spin-wave 
spectrum of amorphous ferromagnets. 

In the present paper we estimate the dependence of structural fluctuations of magnetization on the concentration 
of non-magnetic sites in a dilute ferromagnet. We suggest a self-consistent set of equations for magnetization and its 
quadratic fluctuations derived within the framework of the configurationally averaged two-time temperature Green's 
function method. 



1 Configurationally Averaged Green's Functions 

Let us consider a structurally disordered system of N atoms in the volume V which is described by the isotropic 
Heisenberg Hamiltonian 

H = --^JijmnjSiSj - h^njSj, (1.1) 

i,3 3 

where J,j = J(|R, — Rj|) is the exchange integral describing the interaction between the ith and jtYi atoms, Sj is 
the spin operator of the ith atom, h is the external magnetic field, rii is a random variable taking the value of 1 or 
according to whether or not the site i is occupied by a magnetic atom. 

The Fourier-component of retarding two-time temperature Green's function 

« l\l' » t =« niS+\ni,Sy » t = -»0(f)([S+(t),Sj7(O)]}n,n,, (1.2) 

satisfies within the Tyablikov approximation the equation of motion 

(E-h)< l\l' >b= 26 M /£; + ^2 J iJ ( X J < -xi < j\l' > B ) , (1.3) 

3{&) 

where Xi =< Sf > rii is the magnetization of the ith site and x = ~x~l is the mean magnetization. The overbar means 
the averaging over configurations but due to self-averaging of magnetization we can write 

• T = ^E^ = ^E<^ > ^ = ^ i ^ E <Sf>=n<S?> = na, 

i i i(m = l) 

where double overline means the averaging over the magnetic sites only. Extracting fluctuations = Xj — x and 
coming to a momentum space 



TV , N 

e -iqRj 



« <M » = ± E " H ; « Hj >. e q = -!=5> 

*,i=i i =1 



we get 



(E - E Q (q)) « q|q' »= 2a;5(q - q') + ^=C q - q ' i^E A &> k )^- k « k ^ »' ( L4 ) 

• ' N v N 
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where J(q) = ^ J(fl)e" ,qR is a Fourier transform of exchange integral, Eo(q) = h + xA(q, q) is the spin-wave 

R 

spectrum of a non-dilute crystal and the notation A(q, k) = J(q — k) — J(k) is introduced for convenience. All 
sums in equation ( |l.4| ) and hereafter in the reciprocal space are taken over the Brillouin zone. It is worth while 
remembering that there is just a lattice and therefore all Fourier transformations from the real space to momentum 
one are reversible. 

Note that Xi or £j include fluctuations of structure rij as well as fluctuations of magnetization < Sf >. It is 
convenient to rewrite equation (1.4) in the matrix form: 

G = G° + G°p+G°QG, (1.5) 

where 

Gq.q' q|q' 3> is the unaveraged Green's function, 

G£ „/ = ^ x ^q.q — j s Z ero-approximation Green's function, 
E — E (q) 

p„ „i — ^ q ~ q is the matrix of relative fluctuations of magnetization and the matrix 
q,q VNx 



Qq.q' = — *2 — Pq.q' is proportional to relative fluctuations of magnetization and depends on the exchange 
interaction between spins. The zero-approximation Green's function of the dilute system has the same form as in 
the non-dilute case with the only difference that the mean magnetization x = na must be substituted in place of a. 
Averaging equation (1.5) over configurations we obtain 

G = G° + G°QG. (1.6) 

To obtain an equation for the new unknown Green's function QG we need to multiply equatio n (|l.5| ) by the matrix 
G°Q and to average the product over configurations. Substituting the result in expression we get another 

expression for G: 

G = G° + G°QG°p + G°QG°QG, (1.7) 



which contains the matrix QG°QG. We can obtain the equation for this averaged product in the same way and 
continuing this iteration procedure we can get similarly as it is described in |2j the expression for the configurationally 
averaged Green's function 



G= Jl-G^AQiQ] G° (l + 5^AQip) , (1.8) 

where 



AQi = Qi - Qi, Qi = QG°, Q» = AQj_iQG° for i>l. 
All three factors in equation (|l.8|) are diagonal matrices in the momentum space and we can rewrite this expression as 



G°] q i + e [aq.p; 



[G]q = * (1-9) 

l-[G°] q £[AOQ] q 

i=l 

where [. . .] = [. . .] is a diagonal element of the matrix. 

Let us consider other Green's function Gi = Gp which we will use further. To obtain the equation for this function 
let us multiply the equation (|ll^) by p. We get 

Gx =G°( /9 + p 2 ) + G°QG 1 . (1.10) 

An averaged value of Gi obtained within the similar iteration procedure has the following form 

fl-G fjAQ^ G° (V + fjAQ^p + p^ (l.n) 
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or 

[G°] q ([^] q +E[AQ,(P + P 2 m 

i—1 



[Gi] q = ^ ^5 '- (1-12) 

l-[G°] q £[AQ 4 Q] q 



for the diagonal elements in the momentum space. 

Note that the matrix p 2 is proportional to the unit matrix: 



[p2]q ~ N ^ ^ J q " N S fa - f 2 • (1-13) 

k i 

2 Equations for Magnetization and its Structural Quadratic Fluctuation 

As was shown by Kaneyoshi in Q the averaged Green's function G can be expressed approximately in terms of mean 
magnetization x and its mean quadratic fluctuations £ 2 . To obtain equations describing these quantities let us consider 
the equation for a mean moment at the ith site that expresses it in terms of non-averaged Green's function in the 
energy representation <C l\l ^$>e- In a standard way we obtain 

oo 

Xl = \ n ~ J e^T (~i 3<< 111 >>s +*°) • (2 - 1} 

— OO 

Averaging this equation over configuration we obtain the equation for mean magnetization 

oo 

1 f dE 



X= 2 n - / ^mU 3SpG(£ + !0)| ' (2 - 2) 



where we took into account that the averaged one-site Green's function does not depend on the number of the site: 

« W » = ^ J2 « l \ l » = n H < ^ » = ^ Sp ^- 

I q 



Subtracting the averaged equation (2.2) from the non-averaged one (2.1), multiplying the result by the fluctuation £, 



and averaging the product over configurations we obtain the following equation for quadratic fluctuations 

oo 

T7T 1-. r- f dE 



$ = §An,6 " / ( W »£+«> ) ■ (2-3) 



- oo 



It is easy to prove that 



Ani^i = (ni - n)(ni < S[ > —x) = (1 - n)z 

and 

ii « l\l » = ^SpGT. 

Substituting these two equations in ( |2.3| ) we obtain finally 

oo 

i, s r dE ( i 



;*=x I ~{l-n)- J ^—\-—QSp Gl (E + iO)) | .2.4) 

v — OO 

Let us consider the meaning of mean quadratic fluctuations £ 2 more closely. We find that 

£2 = ( n; < s? >)2 - x 2 = n< Sf > 2 - x 2 = n (< S? > 2 -<T 2 )+ x 2 = n(Acr) 2 + x 2 , (2.5) 
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where the first term in the final expression corresponds just to the fluctuations of mean moment of magnetic atoms 
while the second term corresponds to the fluctuations of structure. Therefore it is convenient to divide these two 
terms. We can easily find that the mean quadratic fluctuations of moments of magnetic atoms have the form 



1 



(Aa)* = a \-(l-n) 



1-n 



dE 



1 



3" E -1 irN 



r 3?SpGi(£; + iO) 



(2.6) 



We can use the equation for mean magnetization (|2.2|) to rewrite this equation in the following form 



oo 



where we have introduced the new Green's function 



C i i — G — C i 



The averaged value of this Green's function can be derived from equations (1.8) and (1.11) as follows 

-l 



GT 2 = (i-G°£aW) G^-^ + i-^f^AQ^-f^AQ; 

\ i=l / \ i=l i=l 



(P + P 2 ) 



or 



(7T — oo oo \ 
+ ^ E [AQ,p] q - E [AQ,(P + P 2 )] q 
i=l i=l / 
L^Jq - 55 



1 - [G°] q E [AQ 4 Q] q 

i=l 

for the diagonal elements in the momentum space. Here we have used the fact that 



-2 f (Aa)^ 1-n 

P = — = — a 

ar no 



ii 



It is convenient to come to the dimensionless variables 

E-h - , s - ft 
£ = ^tt^, /3 = /3J(0), ft 



G Q = J(0)G a . 



(2.7) 



(2.8) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



x J(0) ' " '* J(0)' 

Then the equations for the magnetization x and the mean quadratic fluctuations of magnetization of magnetic sub- 
system (Act) 2 take the form 



x = - | 1 + 2 



d£ 



s f3(£x+h) _ i 



no* 



= 2 



d£ 



e P(Sx+h) _ l 



92(£), 



(2.13) 
(2.14) 



where 



9<*{£) 



2ttN 



3SpG Q (£ + i0). 



(2.15) 



We can assume that the spectral density gi£ ) does not depend on the mean magnetization x, similarly as in the case 
of the non-dilute system. Then substituting zero external field into equation (2.13) for magnetization x and looking 
for the limit of the zero magnetization we obtain in a standard way the following equation for critical temperature 



J(0) 



n / 4 



d£ 



<?(£)■ 



(2.16) 
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To solve this set of equations self-consistently we need to express the averaged Green's functions G and G2 in 

terms of the magnetization x and quadratic fluctuations of the magnetic subsystem only. In the following 

section we consider the limit of low concentration of nonmagnetic impurities satisfying this condition. The analysis 
of different approximations allowing to reach this aim for any concentration of nonmagnetic impurities will be the 
subject of further works. 



3 Low Concentration of Nonmagnetic Impurities 



Let us estimate the value 



1 



Pq,kiPki,k 2 ■ ■ ■ Pk m _i,q 



N m x m 



J2 cx p [ -* E *j ( k J ~ k J-! ) ) 61 • • ■ & 



h,...,l 



(3.1) 



k =k m =q 



we need to calculate to obtain the averaged Green's functions. For the low concentration of non-magnetic atoms we 
can neglect the correlations of fluctuations on different sites supposing 



Hi ■ ■ ■ Hj — Hi - ' ' Hj ■ 

The only values we need to calculate within this approximation are It is easy to show that 



Q n = (1 - n){-x) m + n(Aa + (1 - n)a) m = (1 - n)(-x) m + n{Aa) m + o(l - n). 
Let us remind that for the Gauss distribution of the random value £ 



(AC) 5 



(2m- 1)!!(AC) 5 



(AC) 2m+1 = 0. 



(3.2) 



(3.3) 



(3.4) 



We don't now the distribution law for magnetization fluctuations but we can suppose that the other than quadratic 
fluctuations are small enough to be neglected within the linear as to the concentration of non-magnetic atoms approx- 
imation 

(Act)™ = o ((Act) 2 ) , m > 2. (3.5) 

Hereafter we shall neglect higher than quadratic moments of magnetization of the magnetic subsystem. Within the 
made approximation we obtain 



1 



Pq,kiPki,k 2 • ■ ■ Pk m _!,q 



N 



'71 — 2 / 4 \ j. 



i=2 



x l V x 



3-1 



y~] 6k q - k m _i + ^(k A/ - k A/ -i) 



+ 



k =q 



(3.6) 



where {A}^ n = {Ai , . . . , \j} C {1, . . . , m} and Sk(x) = 6^,0 is Kronecker's symbol. The first term in the right hand 
side of expression (3J3) corresponds to the one-site approximation and the second one corresponds to the two-site 
approximation where we have neglected possible correlations of the fluctuations on different sites. 

Let us calculate the Green's functions G and G2 taking into account only the ter ms linear as to the concentration 
of nonmagnetic atoms. Within the made approximations the series in expressions (1.9) and ( 2.10 ) for the averaged 
Green's functions [G] q and [G2] q can be summed up. Indeed 



£[AQ. iP ] q = £[(QG°)V]q + o(l - n) 

00 

^7V l -^q,ki -B kl , k2 • • ■ -B ki _ liki Pq,kiPki,k 2 • ■ ■ Pki,q 



(3.7) 



EE ! 



i=i 



k i=l 



i+1 



\B' 



q,k? 
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where 

g q ,k = ^ g ^ o k (k) , A q , k = J(q-k)-J(k), £ (q) = 1 - J(q), J(q) = ^j- 
Taking into account the expression for the averaged fluctuations £ m (fTq) we get 



5][Aap] q « ( i-„)(_i)i+i + „L_i_ pv« (3-8) 



k i=l 



(1 - n) £[B(1 + B)^] q , k + £ B qjk . 



k 

In a similar way we can estimate other series as follows 



|>Q 4 Q] q « i ^(1 - n)[B(l + B) _1 A] q . q + ^l![BA] q , q 



(3.9) 



and 

oo 

£[AQ^] q « -(1 - n) £>(1 + B)- 1 ]^. (3.10) 

i=l k 

Thus we obtain the approximat expressions for the averaged Green's functions G and G2 in terms of the mean 
magnetization x and structural magnetic fluctuations of magnetic subsystem only. 



[G] q = 2 k (3.11) 

£ - fo(q) - ^[BA] q , q - (1 - n)[B(l + B)-iA] q , q 



and 



[G 2 ] q = 2 == ^ * i . (3.12) 

£ ~ £o(q) - ^[BA] q , q - (1 - n)[B(l + B)-iA] q . q 



Substituting expression ( ft.!2p in equation ( |2.14| ) we can find that it gives us zero solution for the fluctuations - — V—- 
It means that this value is of a higher order as to concentration of non-magnetic impurities and it should be neglected 
within the linear approximation. Thus we can rewrite the Green's function G in the form 

iei a =2 - l]^-;f {£ :$ £ y (3.13) 

£ - to(q) - (i - "■)£(£> q) 

where 

C(£,q) =^[B(l + B)- 1 ] qik , (3.14) 

k 

S(f,q) = [B(l + B)- 1 A] q , q . (3.15) 

This expression contains neither magnetization nor its fluctuations. The problem of calculating the inverse matrix 
(1 + B) _1 is well known in theory of crystals with impurities and similar matrices were calculated in [0 and ||. We 
present a simple way of calculating the functions C and E for the d-dimensional simple cubic lattice with the nearest 
neighbours interaction in appendix A. 

From the poles of the averaged function G we can obtain the equation for the spectrum of spin excitations 

£-£b(q)-(l-n)E(£,q)=0 (3.16) 



G 



We present here only the long-wave solution of this equation for the simple cubic lattice with the nearest neighbours 
interaction 



E « Dq 2 = xJa 2 1 - 2(1 - n 



I s2 (d) 



aJa'n 1 - 2(1 - n) 



l-I s 2(d) 



q z = aJa z 1 - (1 - n 



1-I s 2(d) 
1 + I s 2(d) 2 



(3.17) 



l-IAd) 



+ o[(l-n) 2 ])q 2 , 



where J is the exchange integral for the nearest sites, a is the lattice constant, d is the dimensionality of lattice. The 
integral 

TT 7T 

1 f . f . sin (ari) 



I s2 (d) 



(2ir) d 



dx„ 



(3.18) 



E (1 - cos(x M )) 



is smaller than the one for all the dimensions except d — 1. In the one-dimensional case the denominator of the second 



term in the expression for spectrum (3.17) diverges that corresponds to the fact that spin- waves can not propagate in 
the one-dimensional ferromagnet with the nearest neig hbours interaction at any dilution. For the dimensions d > 2 
the spin-wave stiffness constant D in expression ( 3.17 ) tends to zero while the concentration of the magnetic sites 
n decreases to the percolation limit n c which within the linear as to the concentration of nonmagnetic impurities 
approximation takes the form 

l-TAd) _ 2IAd) (319) 



Note that we can reduce integral ( 3.18 ) to the one-dimensional form 

^h{t)It\t) 



Is* (d) 



t 



dt. 



(3.20) 



where we have used the transformation 



J2 (1 - cos(x M )) 



and the following properties of the modified Bessel functions l v (i) 



J exp ^-*y)(l - cos(x M ))^di 



(3.21) 



1 

2^ 



: cos(nx)dx — I n (t), 



1 

27 



icosi • 2 , Io(t) - h(t) 

e sin x dx = 



hit) 



(3.22) 



As the expression for the percolation threshold ( 3.19| ) is obtained within the linear over 1 — n approximation it shall 
give reliable results for low dimensions when 1 — n c is small but we can not expect any reliability for large d when 
(1 — n c ) — * 0. Indeed, for d — * oo the percolation threshold n c — > \ that does not coincide with the Bethe expression 
= . 2 d-i w hich is an asymptotic of the percolation threshold at large d. Nevertheless as we can see from the 



Table I the results for the percolation threshold obtained from (3.19) are better for d = 2 and d = 3 than the ones 
obtained within the Bethe-lattice approach. One can find some results of the latter approach for the theory of dilute 
ferromagnets in jl(| O, p"2[ . 

Table I. Percolation thresholds of the sc lattice from this work n c in comparison with the "exact estimates" nj: taken from 
]l3| ] and Bethe expression nf = l/(q — 1). For the sc lattice the number of the nearest neighbours q = 2d. 
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Dimension 


n c 

'"c 




n c 


d = 1 


1 


1 


1 


d = 2 


0.5928 


0.(3) 


0.53306 


d = 3 


0.3116 


0.2 


0.34689 


d = 4 


0.197 


0.(142857) 


0.25585 


d = 5 


0.141 


0.(1) 


0.20286 


d = o 


0.107 


A / AA\ 

0.(09) 


0.16824 


d= 7 


0.089 


0.(076923) 


0.14380 


d = 8 




0.0(6) 


0.12561 


d = 9 




0.058823 


0.11153 


d = 10 




0.052632 


0.10030 


d= 100 




0.0050251 


0.01000 



Let us write an explic it exp ression for the spectral density g(£ ) which we need to calculate magnetization ( 2.1 3| ) 
and critical temperature (2.16) 



(3.23) 



where 



fp n \- 1 cvrri (P x rrfl 1 (1 - n)S"(g,q) + F(g,q) 

o(£, q) = — — ^ tr L(£ + z0) = 5 5-. 

27T 7r(£_£ ( q) _(l- n )S'(£, q )) 2 + ((l-n)E"(£,q)) 2 



F(£, q) = (1 - n)C"(f , q) (£ - f (q) - (1 - q)) + (1 - «) 2 C'(£, q)E"(f , q). 



(3.24) 



The functions C", £' and C", S" are respectively the real and the imaginary parts of the corresponding functions 
( |1J ) and ( ^15|) . They are defined as follows 



C{£ + i0, q) = C(£, q) - , q), £(£ + *0, q) = E'(f , q) - iE"(£, q). 



(3.25) 



We can see that the expression for the obtained spectral density g{£) does not depend on the magnetization and its 
fluctuations. It confirms the assumption we have made in the second section to obtain expression (2.16) for the critical 
temperature. 

Let us estimate now the structural fluctuations of magnetization of a magnetic subsystem at low concentration of 
non-magnetic impurities. To do it we assume that these fluctuations are proportional to the square of concentration 



of non-magnetic atoms and we shall keep in the numerator of the expression ( 2. 1C ) for the Green's function G2 all 
the terms up to the second order as to the concentration (1 — n) of non-magnetic sites. 
Thus we need to sum up the series 



£[AQ,(p + p 2 )] q = £[Q 4 (p + p2 )]q _ £[Q 4 ][p2 ]( 



(3.26) 



i=l 



i=2 



1 - £[(QG°) J ']q £[(QG°)*(p + p 2 )] q - £[Qi][p 2 ] q + o[(l - nf 



J=2 



As we can see from expressions (3.8) and (3.10) the series Yl [(QG°) l (p + p 2 )] q is of the same order as the fluctuations 

i=l 



and therefore it is within the made assumptions of the quadratic order with respect to the concentration of 
non-magnetic impurities. Thus we have 



^[AQ,(p + p 2 )] q « £[(QG°)*(p + p 2 )] q - £[Q,:][p 2 ]c 



(3.27) 
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q.k 



Pq.kPk.q + E Pq.kPk.k'Pk'.q 
k' 



E^ E 5 q,ki-Bki,k 2 • ■ ■ -Bki_i,ki 
i=2 ki,...,kj 



Pq,kiPki,k 2 • ■ ■ Pk;,q + E (Pq,kiPki,k 2 ■ ■ ■ Pk l+1 ,q - 6k;.q Pq.kiPki.ks ■ ■ ■ Pki_i,q Pq,k i+1 Pk i+ i,q) 



n 2 , cc 



E B ^ + E E b im b ) 



J k 



i=2 ki,...,ki 



ki,k 2 ■ • ■ Acj_i,ki x 



7 .7=2 ^ ' ^ / {A} 3-1 \ /=! 



AT ^ 
k i+1 



i+2 



i / £ \ ■? / £ \ ''+ 2 ~J 



J=2 



x / \ x 



3-1 



Mq-k i+ i + ^(k A/ -kx r i) 



i=l k 



i + 2 



+ E E ^q.ki-Bk!,k2 • • -^kj.i.ki X 
i=2 ki,...,k; 



''-I / , \ 3 / £ \ i+l-j 



"Eli 

J=2 



^ 6 q-ki + £(k A ,-M 

{Air 1 V f=1 



k, + i j=2 

It is easy to show that 



EE(Cf (if" E «(q-^ + E(^-^-o) -»(§)' (!)'«*- 



3-1 



E E 6 k ° _ ki +! + E^/ _ kA /- J 



kl + 1 Mill 



/ = 1 



E E 5 ( k o-ki + E( k A/- k A/-i)] +E E 6 (k -ki +1 + ^(k A/ -k A/ _ a 

i-2 



(3.28) 



n J2 6 k -k 4 +^(k A/ -k A/ _o ] 

{AH" 2 V /=1 



where C™ = — rr^ — tt is the binomial coefficient indicating the number of different m-element subsets of the set 

n ml[n—m)\ ° 



{1, . . . , n}. Substituting ( 3.28 ) in the expression ( 3.27 ) we obtain 



(3.29) 



i=i 



_ n ^r + (r>i^g(r(ir^, 



EEt Bi k> 



9 



^E E £ q , kl ...B k< _ liki £ 

i=2 ki ki j=l 



6 q-k l + ^(k A/ -k A/ _ x ) 



k =q 



Now we substitute in this expression the value of y^j from equation (3.3) and finally we obtain 

oo 



+ (l-n) 2 J^5 q , k + (l- 
+ (l-n) 2 )^B q , k + (l- 



where we have used the identity 



2 l - 2. 



Finally we get the following expression for the Green's function G 2 



]r£[B*] q , k (-ir(2 l -2)= (3.30) 

i=2 k 

5>B 2 (1 + 2B)- 1 - 2B 2 (1 + B) _1 ] qi k, 

k 

(3.31) 



[G 2 L = 2- 



^ 1 + £ S q ,k + (1 - nf £[B + (1 + B)- 1 - (1 + 2B)- 1 ],,,, 



£-5>(q) 



(3.32) 



where we have neglected the contribution of disorder to the denominator because of both the terms in the numerator 
are of the order of the square of the concentration of non-magnetic impurities. It is convenient to distincguish the 
two parts of this Green's function as follows 



G 2 



(1 - n) 2 G 2Q - ^=^G 



(Aa) 



2b- 



where 



- C a (£,q) - . 1 + C 6 (£,q) 

l^2aj q - 2- [<J2bJ q = 2- 



£-£b(q)' ^ f-fo(q) 
C„(£, q) = J] [B + (1 + B)- 1 - (1 + 2B)- 1 ] q , k , C b (£, q) = ^ £J q , k . 



(3.33) 

(3.34) 
(3.35) 



k k 

Then we can rewrite the equations for the mean quadratic fluctuations of the magnetization of the magnetic 



subsystem (Act) 2 ( 2.14 ) in the following form 



d£ 



^(Sx+h) _ i 



92a(£) 



1 + 2 



d£ 



e /3(£x+h) _ l 



92b{£) 



where 



92a{£) = _ '^ X 7 SpSG2a ' 92b{£) 



2ttN 



Sp3G 2b . 



(3.36) 



(3.37) 



Let us consider the low-temperature behaviour of structure fluctuations of the magnetization of the magnetic 



subsystem ^ A J$ . Both the integrals in equation(3.36) tend to zero at the temperature tending to zero. Therefore 
we need only the low-energy behaviour of the spectral density <? 2a to obtain the low-temperature asymptotic of the 



fluctuations 



(Act) 2 



In the case of the simple cubic lattice the simple calculations shown in the appendix B yield 



92 



(3.38) 
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giving in the three-dimensional case 



{ ^-^-nfZ mm {^) V \ (3.39) 

where 

oo 

Z v (x) =Y j n- 1/ e~ nx 

n=l 

and J is the exchange integral for the nearest neighbours. 

4 Conclusions 

We have shown that the configurationally averaged Green's function method is a useful tool for deriving self-consistent 
equations describing magnetization and its mean quadratic fluctuations caused by disorder. In particular, at low 
concentration of nonmagnetic impurities we obtain explicit expressions for the spectral densities g{£) and g2{£) giving 

the value of magnetization x and its quadratic structural fluctuations (Ac) 2 with respect to the linear concentration 
of impurities. 

We have revealed that the relative quadratic fluctuations of magnetization ( A J$ could be neglected within the 
linear as to the concentration of nonmagnetic impurities 1 — n. Thay are found to be quadratic over the 1 — n and 
indicate the T 3 / 2 behaviour at low temperatures. 

The same approach can be applied to describe amorphous and liquid spin systems. It will be the subject of further 
works. 



A Calculation of the Functions C and E for the Simple Cubic Lattice 

It is easy to show that the matrix element 



B, 



1 A 



q,k 



q.k 



could be rewritten as 



£ J(R) [(cos(qR) - 1) cos(kR) + sin(qR) sin(kR)] 



B, 



q,k 



R 



N 



£ - 1 - E •W cos(kR; 

R 



For the nearest-neighbour interaction it takes the following form 

. \ E [(cos(ql) - 1) cos(kl) + sin(ql) sin(kl)] 
1 i 



B, 



N 



£- (i- jEcos(ki) 



(A.l) 



where 1 goes over the z nearest neighbours of some site. The fact that the matrix element i? q ,k is invariant under the 
transformation 1 — ► — I allows us to reduce all the sums in expression ( A.l ) as follows 

i i 

where a pri me d enotes that the sum is taken over the z/2 nearest neighbours of some site within one half-space. 
Expression (A.l) can be easily factorized as follows 

BqM= a iM(q)A^(f,k), 



£1=1,2 
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where 



2 2 
ay i(k) = - (cos(kl) - 1) , ay 2 (k) = - sin(kl), 

z ' z 



(A.2) 



A,i(£>k) = 



cos(kl) 



iV 



f - (l-f£cos(kl) 



A, 2 (f,k) 



sin(kl) 



iV 



f- (i-|Ecos(ki; 



Now we can rewrite -B q ,k as a matrix product 



5 q ,k = a T (q)/3(£,k), 



where 



aii ,2 



Ax, 2 



(A.3) 



After the made transformation the inverse matrix (1 + B) _1 can be easily calculated. Indeed 

oo 

[(l + B)- 1 ]^ = 5q , k + ^[(-B)i q , k = 

i=l 

6 q , k ~ « T (q) f 1 + £ (~ E ^ > k ')" T (k')^) J /?(f , k) = 



(A.4) 



where 



6 q , k -a T (q)(l+7(£))- 1 /3(£:,k) ) 



(A.5) 



Thus the problem of calculating the inverse matrix (1 + B) 1 of the order N is reduced to the calculation of inverse 
matrix (1 + 7(f)) -1 of the order z/2. 

Since a± and 0i are even functions of k and ai and p2 are the odd ones the matrix 7 has the form 



m = 



71(f) 
72(f) 



where 



k 

Now we can easily express the functions C defined in ( |3.14| ), ( 3.35|) in terms of the matrices a, (3 and 7 

C(f,q) = E[ B ( 1 + B ) _1 ]^ = a T (q)(l+7(f))" 1 ^(f), 

k 

C„(f,q) = E[ B + ( 1 + B ) _1 -( 1 + 2B ) _1 ]^ = 

k 

d T (q) [1 + 2(1 + 27(f))- 1 - (1 + 7(f))- 1 ] B{£ ), 

c 6 (f, q ) = ^v-« T (qW). 

k 



(A.6) 
(A.7) 

(A.8) 
(A.9) 

(A.10) 
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where 



£(£) = £/3(£,k). 



(A.ll) 



To obtain a similar expression for the function X we rewrite the matrix A in a similar way as the B: 

^ q ,k - a T (q)??(k), 

where 

/ ^1,1 \ 



1 



^7i,i( k ) = cos(kl) T7i, 2 (k) = sin(kl). 



(A.12) 



Thus we can rewrite the expression for the function S as follows 

£(£,q) = [B(l + B) -1 A]q >q = 

« T (q)(l + j(£ ))- 1 7(f ))r)(q) = a T (q)r?(q) - d T (q)(l + 7(^)) _1 ??(q)- 



(A.13) 



The calculation of the inverse matrixe (1 +*/(£)) 1 can be easily carried out for any simple lattice. Let us consider 
as an example the (i-dimensional simple cubic lattice. Then we have 

(kl) = kid, i = 1, . . . , d, 

where a is a lattice constant. The matrix 7 takes the form 

/ <t> X ■■■ X 
X "'■ "'• : 



7 







X 



X 



X 



ip 

'•• 



V 



(A.14) 



... ib / 
All the diagonal elements of the matrix 71 read 

i_ n y cos 2 (fcio) - cos(fcio) _ 1 , d d v 



where as the off diagonal elements read 

_ 1 1 V- cc-sQia) cos(fc 2 a) - cos(fcia) _ 1 , , _ d . 



The matrix 72 = tpl, where 



Here 



£ -£n S 



ib - -- V Sln2(fclQ) - I/ d 



d 

£ sc (k) = l--^cos(M 

i=l 



(A.15) 



(A.16) 



(A.17) 



(A.18) 
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is a dimcnsionless energy of spin excitations of the non-dilute simple cubic ferromagnet and we have introduced the 
following notation 



T d (c) _ 1 V cos 2 (fcia) d _ 1 \ - cos(fcia) 

c2[ ' N L> £ -£g c (k)' cl ' N L> £ -£g c (k)' 

/-d (f) = — V sin2 ( fc i g ) 
s21 J A^^-^W 



(A.19) 



1 ^-v cos(fcia) cos(fc 2 o) 



£-£ s c (k) 



Let us note that the sum ^ (3 that appeared in the expressions for the functions C can be rewritten as 

k 



£/M£,k)=/ c d (£) : , £/3 2 (£,k) = 0. 



(A.20) 



It can be easily verified directly that 

1/d 



(l + 7i)" 



l + 0+(d-l)x 



l/d 



1 + 



/d- 
-1 

V -i 



1 -1 



Now we can rewrite the expressions for the functions C as follows 
(7(5, q) 



-1 \ 



-1 



-1 d-1/ 



l + 0+(d-l)X £-737 



7m 



(A.21) 



(A.22) 



C„(£,q) - af(q)[l + 2(l + 27 1 (£))- 1 -(l+7 1 (£))- 1 ]S 1 (£) = 



£o c (q) 



l + 2^ + 2(d-l)x l + 0+(d-l)x 
1 1 



F - 1 



(A.23) 



C 6 (f, q) = df(q) J] k) - / c d W(q), 



(A.24) 



where we have used the fact that 



1 



1 1 



^ ]v£dX C0S ^ = ^ 

(5-1)^(5) + 3^(5) " ' 



d 



d 



and 



+ (d - 1)X = ^ (*) + V 17 -^ " = 
We can obtain the expression for the function £ in a similar way. 
£(£,q) =a T (q)(l+ 7 (f))- 1 7 (£)) ? 7(q) = 

(q)m(q) - (q)(l + T^))" 1 ^) + 3-^*2 (<0fr(q) = 



(A.25) 



(A.26) 
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1 l d 2(£) l 

-^(cos(fc,;a) - l)cos(fc ia ) + d+ S jd (-gw ^ sin2 ( fc » a ) 



1 1 



l-£I d (£) d^ H( cos ( fc ' a ) ~ !) c ° s (%«) - 

^^(cos(fcia) — 1) cos(fcia) — — ^^(cos(fc;a) — 1) cos(fcj-a) J = 



£b(k)-2ft(k) +arf ^(k) 
(£ (k)-£ 2 (k)) + 



1 - £Ii{£ ) 



where 



(2£ 1 (k)-2£'o(k)+£: o 2 (k)), 



(A.27) 



(A.28) 



Let us consider the long-wave solutions of equation (3.1£) for the spin excitation spectrum. Taking into account 



that 



we can easily find that 



where 



£o(k) = — k 2 + o(k 2 ) and ^(k) = — k 2 + o(k 2 ) 



E(£,q)=-2 



a 2 I_Ad) 
2d 1 - I S 2 (d) 



k 2 +o(k 2 ), 



sin 2 (fcia) 



£\(1 - cos(fcia)) ' 



(A.29) 



(A.30) 



(A.31) 



B Calculation of the Spectral Density g^ a for the Simple Cubic Lattice 

Let us calculate the function 



92a{£) 



2ttN 



Sp3G 2Q 



(B.l) 



we need to obtain the low-temperature behaviour of the structural quadratic fluctuations of magnetization (Act) 2 . 
Using the result of the previous section we get the ci-dimensional (sc) lattice 



92a{£) 



1 So3G~- 1 SV^M 
— Sp^G 2Q - f-f (q) 



2ttN 



^ c (q) 



(B.2) 



£+i0 



i d c {£) 



F - 1 F _ 1 

° 2I%{£) ° Tf{£) 



£+i0 



E-hVe £-v £ 



£+M 



Here 



^ c (q) 



x^- >• ->?; i q) ~^ S) 7 " (£) ' ldl[E) ~ N^£-£ s c (k)> Vs ~W)- 



(B.3) 



15 



Simple calculations give us the following expression for the spectral density 



if(£) + 



(2£ - v^Y + (V^Y {£ - v £ Y + (v^Y 



u'J 



-it (£) 



IV 



V. 



(2£-v>y + (v>'Y {£-vtf + iygy 



where 



Us+io = U' e - iUg, Vs 



my ' £+i0 i*\£)-iif(£) 

Let us consider an asymptotic behaviour of the spectral density g2a(£) at small £ 



V'e + iV'i. 



N ^ £ -£^ c (k) 



£-£g c (k) £^o" 



U 'e = ~Jr E f o C (k)6 {£ - £o c m = -n£g (£), 



V ^ l-£g°(k) 



N^£-£°%k) 



lf{£) = jj E (! - f o c (k)) § (£ £o c (V) = tt(1 - -^So^), 



if(£) 



Here 



(lf(£)) +{lf{£))" (lf(£)Y + (lf(£)) 

•9o( f ) = ^E 5 ( f - f o c (k)) 



2 • 



is the density of states of a non-dilute crystal. Finally we get 



(B.4) 



(B.5) 

(B.6) 
(B.7) 

(B.8) 
(B.9) 

(B.10) 
(B.11) 



\ fa !?(£)+ if (£)U' £ ) - Q 2. 9 o(f). (B.12) 
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